EXACT WAVE FUNCTIONS FOR GENERALIZED HARMONIC 

OSCILLATORS 
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Abstract. We transform the time-dependent Schrodinger equation for the most general variable 
' quadratic Hamiltonians into a standard autonomous form. As a result, the time-evolution of exact 

wave functions of generalized harmonic oscillators is determined in terms of solutions of certain 
Ermakov and Riccati-type systems. In addition, we show that the classical Arnold transformation 
is naturally connected with Ehrenfest's theorem for generalized harmonic oscillators. 



O 
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^ ! 1. Introduction 

I ' 

■ Quantum systems with variable quadratic Hamiltonians are called the generalized harmonic os- 

g . cillators (see [1], [9], [H], [I5], [19], [22], [26], [33], [35], [37], [IQ], [51], [55], [56] and references 

G^l therein). These systems have attracted substantial attention over the years because of their great 
importance in many advanced quantum problems. Examples are coherent states and uncertainty 

^ ■ relations. Berry's phase, quantization of mechanical systems and Hamiltonian cosmology. More 

. applications include, but are not limited to charged particle traps and motion in uniform magnetic 

^ [ fields, molecular spectroscopy and polyatomic molecules in varying external fields, crystals through 

^ ■ which an electron is passing and exciting the oscillator modes, and other mode interactions with 

^ • external fields. Quadratic Hamiltonians have particular applications in quantum electrodynamics 

Q [ because the electromagnetic field can be represented as a set of forced harmonic oscillators [22]. 

y—( A goal of this Letter is to construct exact wave functions for generalized (driven) harmonic 
oscillators [1], [7], [26], [33], [35], [51], [55], in terms of Hermite polynomials by transforming the 
time-dependent Schrodinger equation into an autonomous form [56j. The relationship with certain 
Ermakov and Riccati-type systems, which seem are missing in the available literature in general, 
are investigated. A group theoretical approach to a similar class of partial differential equations 
is discussed in Refs. [T], [B], [12], [21], [12], [IS] (see also [IS], [I9j and references therein). Some 
applications to the nonlinear Schrodinger equation can also be found in Refs. 
, [17] and 



2. Transforming Generalized Harmonic Oscillators into Autonomous Form 
We consider the one-dimensional time-dependent Schrodinger equation 

.f^.., (2.1) 
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where the variable Hamiltonian H = Q {p, x) is an arbitrary quadratic of two operators p = —id/dx 
and X, namely, 

= -a (t) + b (t) x^^P - %c (t) xtP, -tdit)^-f (t) xij + ig (t) (2.2) 

(a, 6, c, (i, / and g are suitable real- valued functions of time only). We shall refer to these quantum 
systems as the generalized (driven) harmonic oscillators. Some examples, a general approach and 
known elementary solutions can be found in Refs. [7], [S], 0, HO], [11], [20], EU, [22], [S5], [3S], 
[3S], [IS], [M] and ^55j. In addition, a case related to Airy functions is discussed in [32] and Ref. [TT] 
deals with another special case of transcendental solutions. 

The following is our first result. 

Lemma 1. The substitution 

i(a{t)x^+5{t)x+K(t)) 

^ = r=TT ^(^'^)' ^ = Ht)x + e{t), T = ^{t) (2.3 

V/i it) 

transforms the non- autonomous and inhomogeneous Schrddinger equation Ii2.2\) into the autonomous 
form 

- ^Xr = + CoC'x (Co = 0, 1) (2.4 

provided that 

^ + b + 2ca + Aaa^ = coa/?^, (2.5 
dt 

^ + (c + 4aa)/3 = 0, (2.6 



dt 



^ + a/32 = (2.7 
dt 



and 



Here 



— + (c + Aaa) 6 = f + 2ga + 2coa/3^e, (2.8 

(JjL 

ds 

- = (^7 - 2a6) /3, (2.9 

dK 

^ = gS - a6^ + coaph^ (2.10 
dt 

1 11' d 

a = —— . (2.11 

4a /i 2a 



Proof Differentiating ^ = /x-^/^ gi5(x,j)^ (^^^ ^i^j^ S = a (t) x^ + 6 (t) x + k (t) , ^ = /3 (t) x + e {t 
and r = 7 (t) yields 

1 



{a'x^ + S'x + k')x + i( W'x + e') + YXr - 



(2.12 



e-^'i^^ = ^ (2ax + 5) X + (2-13 

and 

e-*^7/>,, = [(2«a - (2ax + df) x + 2^ (2ax + 6) /3x^ + . (2.14 
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Substituting into 

iijj^ = —ail^xx + {b — Coals'^) x'^ip — icxip^ — idijj (2-15) 

and using system f l2.5p - fl2.10p . results in Eq. (12. 4p . Further computational details are left to the 
reader. □ 

Our transformation f l2.3p provides a new interpretation to system f l2.5p - fl2.10p originally derived 
in Ref. [7] when cq = by integrating the corresponding Schrodinger equation via the Green function 
method (see also [SD] for an eigenfunction expansion). Here, we discuss the case Cq 7^ as its natural 
extension. 

The substitution (12. lip , which has been already used in [7], appears here from a new "transfor- 
mation perspective" . It now reduces the inhomogeneous equation (12.50 to the second order ordinary 
differential equation 

/i" - r (t) /i' + 4(T (t) /i = Co {2af /^V, (2-16) 
that has the familiar time-varying coefficients 

T{t) = 2c + Ad, a{t) = ab- cd + d^ + - { . (2.17) 

a 2 \a d J 

(The reader should be convinced that this derivation is rather straightforward.) 

When Co = 0, equation (12.51) is called the Riccati nonlinear differential equation [S2], [S3]; conse- 
quently, the system fl2.5p -f l^TTU]) shall be referred to as a Riccati-type system. (Similar terminology 
is used in [19] for the corresponding parabohc equation.) Now if co = 1, equation f l2.16p can be 
reduced to a generalized version of the Ermakov nonlinear differential equation (15. ip (see, for ex- 
ample, [9], [18], [M], [50] and references therein regarding Ermakov's equation) and we shall refer 
to the corresponding system fl2.5p - fl2.10p with co 7^ as an Ermakov-type system. 



3. Green's Function and Wavefunctions 

Two particular solutions of the time-dependent Schrodinger equation (12. 2p are useful in physical 
applications. Using standard oscillator wave functions for equation (12. 4p when Cq = 1 (for example, 
[23] , [3T] and/or [?T]) results in the solution 

i(^ax^+5x+iij+i{2n+l)'Y 

(x, t) = 7=^= e-(^^'+^)^/2 {f3x + e) , (3.1) 

where Hn{x) are the Hermite polynomials [13], provided that the solution of the Ermakov-type 
system f l2.5p -f j2rT0l) is available. 

The Green function of generalized harmonic oscillators has been constructed in the following 
fashion in Ref. [7|: 

G (z, y, t) = , ^ , . exp [i {ao (t) x^ + /3o {t) xy + (t) y^ + 60 (t) x + eo{t)y + kq (t))] . (3.2) 
V 27rz/io {t) 
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The time-dependent coefficients a^, (3q, 7q, 6o, Eq, Kq satisfy the Riccati-type system fl2.5p - fl2.10p 
(co = 0) and are given as follows [7], 



«o(t) = y^#l-^, (3.3) 

" ^ ^ 4a (t) ^0 (t) 2a (t) ' ^ ^ 

/3o (t) = A (t) = exp (- f (c (.) - 2d {s)) ds\ , (3.4) 

^°^'^-2/.,(0);zo(t) + 2a(0) ^^'^^ 



and 



/io (t) 



L 



r / \ (-5) / \\ / \ 9{s) , , . 



a (s) I 2a is) 



ds 



(3.6) 



2a (t) A (t) , , , ^ ais) a is) \ {s) , . x , . xx , . „x 

^o(t) = V-f^(5o (t) + 8 / \\ ^ /io (^) ^0 (3)) ds (3.7) 

A^o (/io (s)) 



«0 (^) 



(^0 (») «o is)? is (3.8) 

/^O Jo (/io (-5)) 

-2 4^ (/io (^) 5o (^)) (f is) - is)) ds 

(5o (0) = — £o (0) = g (0) / (2a (0)) and (0) = 0) provided that /Xg and /x^ are standard solutions 
of equation f l2.16p with cq = corresponding to the initial conditions /Xq (0) = 0, /Xg (0) = 2a (0) 7^ 
and (0) 7^ 0, /i^ (0) = 0. (Proofs of these facts are outlined in Refs. [7], [11] and [16]. See also 
important previous works [15], [37], [5l], [55], [56] and references therein for more details.) 

Hence, the corresponding Cauchy initial value problem can be solved (formally) by the superpo- 
sition principle: 

poo 

^{x,t)= / G{x,y,t)ijiy,0) dy (3.9) 



for some suitable initial data ip (x, 0) = ip (x) (see Refs. [7], [16] and [50] for further details). 
In particular, using the wave functions (13. ip we get the integral 

POO 

^Ij^ (x, t)= G (x, t) ^„ (y, 0) dy, (3.10) 



and this can be evaluated by 

/oo 
e-^'i--yfH^{ay) dy (3.11) 
-00 



^ 2\"/2 / AaX I ,2 



which is an integral transform equivalent to Eq. (30) on page 195 of Vol. 2 of Ref. [16] (the Gauss 
transform of Hermite polynomials), or Eq. (17) on page 290 of Vol. 2 of Ref. [17] . 
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4. Solution to Ermakov-type System 



As shown in the previous section, the time evolution of the wave functions (13.11) is determined 
in terms of the solution to the initial value problem for the Ermakov-type system. In this section, 
formulas f l3.ip -( l3l2|) and f l3.10p -( l3.1ip shall be used in order to solve the general system (12. 50 - 02. lOp 
when Co 7^ along with the uniqueness property of the Cauchy initial value problem. At this point, 
we must remind the reader how to handle the special case Cq = considered in |46j . 

Lemma 2. The solution of the Riccati-type system l{2.5\) - l[KT^) (cq = 0) is given by 

/i(t) = 2/i(0)/^oW(a(0) + 7oW), (4.1) 

a (t) = ao (t) - (4.2) 
^' 4(a(0)+7oW)' ^ ' 

/5(t)^-^4^^-^MA(t), (4.3) 
2(a(0)+7oW) /^W ^ ^ 

and 

A m A m /3o jt) (^(o) + eo(t)) 

^^'^ = 2(a(0)+7o(t)) ' ^'-'^ 

e{t) = £ -r— , 4.6 

2 [a (0) +7o (t)) 

in terms of the fundamental solution l{3.3\) - l[37^) subject to the arbitrary initial data /U (0) , a (0) , 
^(0)^0, 7(0), 5(0), ^(0),«:(0). 

This solution can be verified by a direct substitution and/or by an integral evaluation. This result 
can also be thought of as a nonlinear superposition principle for the Riccati-type system and the 
continuity with respect to initial data holds ^46j . 

Hence, the solution (14. II) - (14. 71) implies the following asymptotics established in ^B] : 

Ko(«) = 0(i) (4.12) 
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as t — )■ for sufficiently smooth coefficients of the original Schrodinger equation f l2.2p . Therefore, 



G{x,y,t) 



1 



^27Tia (0) t 



exp 



' 4a (0) t 



(4.13) 



X exp 



—I 



a' (0) 



{x - y f + 



c(0) 



X 



2\ ^(0) . ^ 

y ) - TTTTT^ - y) 



^a2(0) ' 4a (0) ^ ^ ^ 2a (0) 

as t — > (where / ~ (7 as t — )■ 0, if linit^o (//fi') = !)• This corrects an errata in Ref. [7] . 

Finally, we present the extension to a general case when cq 7^ 0. Our main result is the following. 
Lemma 3. The solution of the Ermakov-type system Ii2.5\) - I[KT^) when cq = 1 (7^ 0) is given by 



/X = /i (0) /io^/?' (0) + 4(a(0) + 7o)' 
«(0) + 7o 
(0) + 4 (a (0) + 7o)" 
/3 (0) /3o 



a = ao — I3q 



(3 



/3(0)/i(0) 



(0)+4(a (0)+7of 



, , 1 /3^(0) 

7 = 7(0) arctan — ; — — -. 

^ 2 2(a(0)+7o)' 



/i(t) 
a (0) > 



A(t) 



6(0) (3' (0) + 2(a(0)+7o) (^ (0) + £0) 



(^0 - /^O 

2£ (0) (a (0)+7o)-/3(0) (5 (0) + £0) 



(0) + 4 (a (0) + 7o)' 



^/3^(0) + 4(a(0)+7o)' 

,3 (0) + ^0 



K (0) + Ko - £ (0) r (0) 



(0) + 4 (a (0) + 7or 



(4.14) 
(4.15) 

(4.16) 
(4.17) 

(4.18) 
(4.19) 

(4.20) 



+ (« (0)+7o) 



(0) (0) - {6 (0) + £0) 
(0) + 4 (a (0) + 7o)' 



m terms of the fundamental solution l{3.3\) - l[37^) subject to the arbitrary initial data /i (0) , a (0) , 
/3(0)7^0, 7(0),5(0),£(0),/t(0). 

Following are the steps to the sketch of the proof. Evaluate the integral fl3.10p with the help of 
f l3.1ip by completing the square and simplify. Use the uniqueness property of the Cauchy initial 
value problem. One can also verify our solution by a direct substitution into the system fl2.5l) - fl2.10l) 
when Co = 1. These elementary but rather tedious calculations are left to the reader (the use of a 
computer algebra system is helpful at certain steps). 

Furthermore, the asymptotics f l4.8p - fl4.12p together with our formulas f l4.14p -f l4.20p result in the 
continuity with respect to initial data: 

lim u (t) = u (0) , lim a (t) = a (0) , etc. (4.21) 
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Thus the transformation property fl4.14p - fl4.20p allows us to find a solution of the initial value 
problem in terms of the fundamental solution fl3.3p - fl3.8p and it may be referred to as a nonlinear 
superposition principle for the Ermakov-type system. 

5. Solution of the Ermakov-type Equation 

Starting from fl2.16p - fl2.17p when Cq = 1, and using fl4. 16p we arrive at 

fi" - r (t) /i' + 4a (t) fi = (2a)' (/3 (0) /i (0) A)^ (5.1) 

which is a familiar Ermakov-type equation (see [S]) 13) [IB], ISO], |SS] and references therein). 
Then our formula fl4.14p leads to the representation 

^ ^ ' - (0) f.lit)+(^ + llMi^) ' (5.2) 



/ a' \ a' 

+ [Aab-c^ + c c' ]x = 2af - g' + g eg. (6.2) 

\ a J a 



^/i(0)y ^ v^^uvy V/ujO) 2/i(0)a(0) 

given in terms of standard solutions /ig and fii of the linear characteristic equation fl2.16p when 
Co = 0. Further details on this Pinney-type solution and the corresponding Ermakov-type invariant 
are left to the reader (see also [5j and [50J). 

6. Ehrenfest Theorem Transformations 
By introducing an expectation value of the coordinate operator in the following form 

one can derive Ehrenfest's theorem for the generalized (driven) harmonic oscillators (see, for exam- 
ple, [8] and [9]). Then the following classical equation of motion of the parametric driven oscillator 
holds 

(fx a' dx f . , o a' ,\ _ „ „ , a' 
dt'^ a dt 

The transformation of the expectation values 

^ = /3x + £, e=(x,ex) with (x,x) = l, (6.3) 

corresponding to our Lemma 1, converts (16.20 into the simplest equation of motion of the free 
particle and/or harmonic oscillator: 

-A + 4co^ = (Co = 0,1). (6.4) 
(This can be verified by a direct calculation.) 

Remark 1. An exact transformation of a linear second-order differential equation into the equation 
of motion of free particle was discussed by Arnold [3J. An extension of the later to the case of the 
time- dependent Schrddinger equation had been considered, for example, in Ref. [56] and recently it 
has been reproduced as the quantum Arnold transformation in [2J and [25J (see also [I], [H], [12], [21]; 
[30] . [12], [15] and [51] for similar transformations of nonlinear Schrddinger and other equations of 
mathematical physics). We elaborate on a relation of the quantum Arnold transformation for the 
generalized (driven) harmonic oscillators with a Riccati-type system when Co = (transformation 
to the free particle) and consider an extension of this transformation (in terms of solutions of the 
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corresponding Ermakov-type system) to the case Cq 
oscillator ). 



= 1 (transformation to the classical harmonic 



7. Conclusion 

In this Letter, we have determined the time evolution of the wave functions of generahzed (driven) 
harmonic oscillators (13. ip . known for their great importance in many advanced quantum problems 
[22], in terms of the solution to the Ermakov-type system (l2.5p -f l2TT0|) by means of a variant of the 
nonlinear superposition principle (I4.14p - (l4.20p . Moreover, the classical Arnold transformation is 
related to Ehrenfest's theorem. Numerous examples, the corresponding coherent states, dynamic 
invariants, eigenfunction expansions and transition amplitudes [U], [52], [S3], [SS], [SZ], [SB], [SO] 
will be discussed elsewhere. 
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